Analysis of exchange terms in a projected ERPA Theory applied 

to the quasi- elastic (e, e') reaction. 



E. Bauer * 

Departamento de Fisica, Facultad de Ciencias Exactas, 
Universidad Nacional de La Plata, La Plata, 1900, Argentina. 



A. Polls and A. Ramos 
Departament d'Estructura i Constituents de la Materia, 
Universitat de Barcelona, Diagonal 647, E-08028, Spain. 
(February 9, 2008) 



Abstract 



A systematic study of the influence of exchange terms in the longitudinal 
and transverse nuclear response to quasi-elastic (e, e') reactions is presented. 
The study is performed within the framework of the extended random phase 
approximation (ERPA), which in conjuction with a projection method permits 
a separation of various contributions tied to different physical processes. The 
calculations are performed in nuclear matter up to second order in the residual 
interaction for which we take a (7r + p)-model with the addition of the Landau- 
Migdal (/-parameter. Exchange terms are found to be important only for the 
RPA-type contributions around the quasielastic peak. 
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I. INTRODUCTION 



The nuclear response to an electromagnetic probe is a common tool to investigate the be- 
haviour of the atomic nucleus [0J . At variance with a hadronic probe it allows a perturbative 
treatment in the external operator coupling constant. In this work we will concentrate on 
the study of the nuclear response function for longitudinal and transverse inclusive quasi- 
elastic electron scattering reactions. These responses are experimentally separated |2| [7| . 
showing that a simple model like the Fermi gas fails to reproduce the experimental data. 
The attempts to go beyond this model can be classified in two groups. On one side there 
are the methods that assume the nucleus as an assembly of non interacting nucleons with 
individual properties, such as the charge radius, modified with respect to the vacuum due to 
the presence of the other nucleons P-P]. Another option is to explore the possibilities of a 
rigorous many-body theory ||l"2|-|36|l keeping the nucleons as the essential degrees of freedom 
with the same properties as in the vacuum, before resorting to such exotic effects. This work 
falls in this second strategy. 

Several approaches to the nuclear many body problem of the nucleus for these processes 
have been extensively analyzed in the literature Microscopic many-body theories must 
deal with short range correlations (SRC) originated from the short range repulsion of the 
nucleon nucleon (NN) interaction. Variational calculations account for SRC by introducing 
a Jastrow correlation factor explicitly in the wave function. In this way, it is possible to 
define a correlated basis function (CBF) and build a fast converging perturbation theory 
using this basis. There are recent studies of both the longitudinal and the transversal 
response for nuclear matter in this framework Hl6| . Alternatively, the effect of SRC can 



be incorporated by introducing a well behaved effective interaction of G-matrix type or the 
standard Landau-Migdal parameterization, with which one can perform perturbation theory 
to build other correlations, for instance of RPA type. As it is not our aim to describe the 
full set of approaches, we will comment only on the ones which lead to our particular theory, 
which is based on the perturbative approach. 

A simple way to introduce in the response the nucleon-nucleon correlations originated 
by the residual interaction is by means of the RPA theory, where one particle-one hole 
excitations are summed up to infinite order. Although improving the Fermi gas picture, the 
RPA approximation is not able to explain some features of the response as, for instance, the 
strength in the "dip" region of the transverse response. 

An improvement upon the RPA-theory consists of allowing the coupling of one particle- 
one hole states to two particle-two hole ones. This corresponds to what has been called final 
state correlations. Two formalisms study this kind of processes. One is the second RPA 
(SRPA) theory (see Ref. [ffH] and references therein) and the other is the Green function 



scheme of Ref. |L7j (see also fi8H2T|). The first one introduces final state correlations over 
the particle-hole bubbles of the RPA theory. In the second one the relationship between 
forward virtual Compton scattering and inclusive electron scattering is used to construct 
a one-body approximation to quasi-elastic electron scattering. In fact, at large momentum 
transfers, where the effect of long range correlations is negligible, the SRPA and the optical 



model Green's function approach should coincide [jig] . 

Both the SRPA and the optical-model Green function approaches use the full residual 
interaction and allow for many particle-many hole final states. Still in both approaches the 
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one body external operator is limited to create (or destroy) a one particle-one hole pair. Once 
the external operator is allowed to scatter a particle (or hole), then two particle-two hole 
states stemming from ground state correlations (GSC) could be activated. The importance 



of these GSC are particularly relevant in the dip region for the transverse channel p3 |. A 



theory to calculate the response function which takes into account all the above requirements 
is already established as the Extended RPA (ERPA) theory [p5| , p7|| . Still, the application 
of the ERPA theory is in general a prohibitively large task. In Ref. |^U| we developed a 
projection method which extracts the main ingredients of the ERPA theory. In that work, 
the scheme was presented and the response was calculated neglecting the exchange part in 
the matrix elements of the nuclear particle-hole interaction. Therefore it seems necessary to 
complete the scheme by investigating and stablishing the influence of the exchange graphs 
in the longitudinal and transverse response of non- relativist ic nuclear matter. Actually, 
the importance of the exchange terms in the RPA theory is a well known problem which 
can not be satisfactorily solved for finite range interactions, although several attempts have 



been made. In a previous work pi| , we developed a simple scheme to evaluate the full 
antisymmetric RPA series contributing to the nuclear matter response, with the result that 
exchange contributions are important, specially at low momentum transfer, and can not be 
accounted for by simply evaluating the ring (direct) series with an effective g' parameter. 

In view of the importance of the RPA exchange terms and the ongoing efforts to improve 
the description of the nuclear response it seems appropriate to explore whether the exchange 
terms in the remaining types of diagrams contained in the ERPA scheme are also important. 
In this work, we undertake this task and evaluate for the first time the contribution of the 
exchange terms to the nuclear response up to second order in the effective interaction. 

The paper is organized as follows. In Sect. 2 we present the formalism. In Sect. 3 the 
results for the exchange contributions to the nuclear matter structure function are presented 
and compared with the corresponding direct ones. Finally, the conclusions are summarized 
in Sect. 4. 



II. FORMALISM 

The Longitudinal (L) and Transverse (T) structure functions, Sl,t{^ hu), are defined 

as 

S L>T (q, M = — Im < 0\O LT ^G(hu)O LT \0 > , (1) 

7T 

where hu represents the excitation energy and q the magnitude of the three momentum 
transfer. The nuclear ground state is denoted by |0 > while the polarization propagator 
G{fiuj) is given by 

G W = n WJ- ~ * A ■ ' ( 2 ) 

nuj — H + it] noo + H — irj 

where H is the nuclear Hamiltonian. Explicit forms for the external excitation operators 
Ol,t are given by 

L = ^l + nU) e iq-x i} (3) 

3=1 
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Ot = ^£{I±|^[g x {p . e ^ ]+i ^±Jp^l [q x (<t(j) x g)] e «^}, (4) 

where m is the nucleonic mass, a^- and p., denote the position and momentum operators for 
individual nucleons, and fi s = 0.88, fi v = 4.70 are related to the proton and neutron magnetic 
moments. In fact, Ol is the charge density operator while Ot is related to the convection 
and magnetization current density. The structure functions are related to the response 
functions, Rl,t(q, hu), through the usual dipole electromagnetic form factor Tiuj): 

We introduce now the usual projection operator P, which projects into npnh configura- 
tions with n — 0, 1, defined respect to the HF vacuum, which corresponds to the case n = 
and is denoted by | >. In addition, two projection operators Q and R are introduced. The 
action of Q (R) is to project onto the npnh space with n being an even integer greater or 
equal to two (n odd greater or equal to three). Explicit expressions are given by 

p= |n><n|, (6) 

n = 0,1 



Q= J2 \n><n\, (7) 
n even 
n > 2 



R= J2 \n><n\ , (8) 
n odd 
n > 3 

where \n > indicates a npnh configuration. 

In the literature only one projection operator, which is the sum of Q and R, is usually 
used. The present separation is done for convenience as it helps to clarify the role of 3p3h 
configurations (see Ref. 0). It is easy to verify that P + Q + R = 1, P 2 = P, Q 2 = Q, 



R 2 = R and PQ = QP = PR = RP = QR = RQ = 0. 
Inserting the identity in Eq. ([I]) one obtains 

S = Spp + SpQ + Sqp + Sqq + Spp t + Srp + Srr + Sqr + Srq , (9) 

where, for simplicity, we have omitted the subscripts L, T. The expression for Spp is given 
by 

/ 777 

S PP (q, hu) = — < 0\O j PG(hcu)PO\0 > ; (10) 

and similar expressions can be written for Spq, etc. To evaluate the propagators PGP, 
PGQ, etc., one has to solve the following equation: 
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Q-Q- 



(11) 



where 



(PGP PGQ PGR 
QGP QGQ QGR 
\RGP RGQ RGR 



POO 
Q 
R 



This is an easy task once the properties of the projection operators are employed. Keeping 
terms up to second order in the nuclear interaction, all 3p3h final state contributions cancel 
each other |30|]. Thus, from all terms of Eq. @, only Spp, Sqp (Spq) an d Sqq survive: 



•>pp 



Im 



71 



hoo-H- - ReY, PRP + i V 



PO\0 >. 



(12) 



QP 



<0O f P 

7T 



htu — H + irj 



P P^resQ 



and 



QQ 



Im 



7T 



< 0\O ] Qr 



hcu — H + irj 



QO\0 >, 



QO\0 >, 



(13) 



hcu — H + irj 

where Sqp equals Spq. The self energy operators introduced in Eq. 

1 



(14) 



are given by 



S — PH res Q 



flLJ — H + 17] 



QH res P 



(15) 



and 



ReZ 



PRP 



V 

P H res R— — — RH res P , 
hoo — H 



(16) 



where V denotes the principal value. We have separated the total Hamiltonian H into a 
one-body part, Hq, and a residual interaction, H res . 

As pointed out in Ref . |3(| , there is still a contribution stemming from a 3p3h configura- 
tion given by the real part of 3p3h- self-energy insertion (ReT, PRP ) in Eq. ([12]). That is, up 
to second order, no 3p3h physical state is possible, but virtual intermediate 3p3h configura- 
tions produce a shift in the ground state energy (see Ref. and Ref. ]2B| for more details). 
The next step is to establish the structure of the ground state. Including the ground state 
correlations perturbatively one gets, up to first order in the residual interaction, 



|0 >= | > -Ho^QHresPl > , 

with | > being the Hartree Fock ground state. More explicitly, one can write 



10 >= 



> 



E 



£ 



|2>, 



(17) 



(18) 



gsc 



where the quantity e gsc refers to the energy of the first order correction to the ground state 
energy. 
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The aim of this section is to present the formalism showing explicitly antisymmetric 



matrix elements. The guidelines to obtain analytical expressions are given in Ref. ||30|| , where 
direct contributions were studied. To complete the scheme, expressions for the exchange 
self-energy insertions are given in Appendix A and exchange terms to the structure function 
are presented in Appendix B. Also, and in order to simplify the calculation, we will limit 
ourselves to the case in which the external operator is attached to the same bubble. We 
will study some exceptions to this as a consequence of antisymmetrization. The three non 
vanishing contributions, Spp, Sqp and Sqq, will be analyzed separately below. Special 
attention will be paid to Spp as its structure is very rich and represents the main contribution 
to the response function. 



A. Spp contribution: 

Let us carefully analyze all graphs stemming from Spp. To do this, we first insert the 
definition of P given by Eq. (||) into Spp [Eq. (|12|)]: 

Spp = Im V < 0\O j \n >< n\ 1 

n,n'=l. 



hu-H Q - H res - S p Q p - ReJ: PRP + irj 



\ri >< n'\O\0 > . (19) 

Using the ground state given by Eq. flTS| ) in the expression and neglecting all third and 
higher orders terms except the ones with self energy insertions, one can write 

g g Lindhard+self —energy _|_ g first order RP A _|_ <j secondorder RPA (20) 

where 

n Lindhard+self -energy _ _^ j . I ( n\ l-i . , -i I \ 

^ PP ~~n Im ^<\ U |1><1| ^_iJ -S^- J ReE^ + ^ 

|1' > a < l'\0\ > , (21) 

Sppfl rst order RPA = _ 1 j m jy < l^th > \ <\\E reb \V > a 

ix Y\, nuj — e\ + ii] 

< l'lOl > 



TibJ — By + it] 



_ 2 yy <\H res \2> a <2 | 0t |l> \ <1|0|>}, (22) 



Spp seeondorderRPA = __ jc y < i^fh > _ < \\H res \V > a 

ix 1 ^ v , Tiuj-ex + irj 
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<C 1 I H res 1 1 ^> a 



Hu> — S\i + it] 
1 

TOO — S\ii + Z77 



- 2 EE <|gres|2>a <2|o f |i> 



< l\H res \l' > a l —— < l'\0\ > 



hco — e\ + irj ' huo — E\> + ir\ 

+EE < lHresl2>a <2|ot|i> 

1 2,2' £ 9SC 



I < 1[0 | 2 /> <2 \ H res\ >a } ? (23) 



fiu-ex+ir) e gsc 

and i?o|w >= £ n |™ >• 

From Eqs. ( p!5|) and the self-energy insertions read now 

< 1|S PQP |1' >= E < MHres\2 > a 1 — < 2\H res \l' > a , (24) 

2 nu — e 2 + it] 

and 

< l\ReZ PRP \l' >=-E< l|gre,|3 > a - - < 3|# res |l' > a . (25) 

3 hcu — e 3 

In all expressions we have made explicit indication of antisymmetric matrix elements. 

In Fig. 1 some graphs contributing to Spp are shown. Within brackets we have col- 
lected direct plus exchange contributions. Let us start by analyzing the contribution to 
g pp L ™dhard+seif -energy _ rpj^ p resence Q f se if ener gy operators makes the energy denomina- 



tors in the r.h.s. of Eq. pi]) to be non diagonal in our particle-hole basis. Non diagonal 
terms, shown by graphs SE2D, SE2E and SE2E' in Fig. 1, are evaluated at second order. 
For diagonal ones, shown by graphs SE1D, SE1E, SE3D and SE3E in Fig. 1, we first 
build up an antisymmetric self-energy insertion and then sum it up to infinite order. 

The first two orders leading to the RPA response (eqs. fl22| ) and (p3|)), are shown by 
graphs RPA1D to RPA2E' in Fig. 1, where only the forward going contributions are 
explicitly shown, that is, the ones stemming from the first terms on the r.h.s. of the above 
mentioned equations. If exchange terms were neglected one would be able to sum terms up 
to infinite order leading to the usual ring series. 

As mentioned, in this paper we keep terms up to second order in the evaluation of 
exchange contributions. However, for the RPA type diagrams we use the method described 



in Ref. which allows to effectively sum up the full antisymmetric RPA series. That 



method is based on splitting the interaction into a pure contact part and a remaining part 
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chosen such that the second-order ring coincides with the full ring series. The pure contact 
term allows a straightforward evaluation of the antisymmetric RPA series up to infinite order, 
while only terms up to second order are retained for the remaining part of the interaction. 



B. Sqp contribution: 

Using the definition of P and Q into Sqp [Eq. (|T3"D1 we have 

S QP = -2- Im{ £ < \0'\l > 1 — < l\H res \2 > a 

71 f^ 2 , TlOJ — E\+ 17] 

I < 2 | | 2 , > <V\Hres\ >a } ^ 



huo — 62 + irj e 



gsc 



Note that as O is a one body operator, it can scatter a particle (or hole) or create (or 
destroy) a particle-hole pair. That is, the Hartree-Fock ground state is not connected to a 
2p2h configuration through O. 

In Fig. 2 we present the second order contributions to Sqp. As a consequence of 
antisymmetrization, a direct term where the external operator is attached to a different 
bubble has come into play (given by graph SQPD' of this figure). Naturally, when we act 
with the antisymmetrization operator over this term, the same set of graphs appears. We 
account for this through a factor two. 



C. Sqq contribution: 

Finally, the expression for Sqq [Eq. (JH|)] is simply 

Sqq = -- Im{ £ < lHresl2 >a < 2|0t| 2 ' > l - < 2'|0|2" > 

7i 2 ^r 2 „ e gsc hu - e 2 > + irj 

< 2"\H re8 \ > a ^ ^ ^ 

£gsc 

Here, the only possible action for the external operator is to scatter a particle or a hole. 

In Fig. 3, we present the main contributions to Sqq. Graphs SQQ3D' presents a direct 
contribution with the external operator attached to a different bubble, in complete analogy 
with Sqp. 

III. RESULTS FOR NUCLEAR MATTER 

In order to benefit from the advantage of translational invariance, results have been 
obtained for infinite nuclear matter at normal saturation density corresponding to a Fermi 
momentum kp = 1.36 fm _1 . 



S 



For the residual interaction H res we assume the ( tc + p )-exchange model at the static 
limit with the addition of the Landau Migdal (/-parameter. In pionic units it reads 



H rea (l) = ^Tl(l)(g'(l)r • r' a ■ a' + ~h\l)r -t'ct-T*'- T) 



(28) 



with 



9[) ~ 9 ~m) p ~p+vy 



(29) 



/ 2 r 2 ^) / 2 

where /x^c (/J, p hc ) is the pion (rho) rest mass and C p = 2.18. For the form factor of the 
nNN (pNN) vertex we have taken 

r (i) - A ^ - { ^ phcf (3D 

* A) ~ A 2 ip + (^c/) 2 ' {6l) 

with A n = 1.3 GeV and A p = 2. GeV. The role of the (/-parameter is to account for short 
range correlations. Note that for a pure contact interaction exchange contributions have been 
traditionally included in the RPA series by a redefinition of the Landau-Migdal parameters. 
In particular, standard (/-values range from 0.7 to 0.95 but, when redefined to account for 



antisymmetric terms, the values are lowered and range from 0.5 to 0.7 [37j. As we evaluate 
explicitly exchange graphs a standard (/-value from 0.7 to 0.95 should be used. We have 
employed the value g' = 0.7. In addition, in all the diagrams considered in our calculations, 
the nucleon lines have been dressed in an average way by taking a momentum independent 
effective mass value of m* jm = 0.85. 

From Eqs. (21)-(]27f)) explicit expressions for the structure functions in nuclear matter 
can be obtained, where the sums over the different configurations are replaced by multidi- 
mensional integrals. The expressions of the direct terms were reported in Ref. |3(| and will 



not be repeated here. In Appendix A we give the exchange self-energy insertions appearing 
in Eq. (|21| ) (see Eqs. (|24|) and (p5|)), while in Appendix B the exchange contributions to 
Eqs. (PTj)-(E3) and (G5p-(E7|) are shown (see also Figs. 1-3). The multiple integrations have 
been performed using a Monte Carlo technique. 

Let us analyse the three non- vanishing contributions, Spp, Sqp and Sqq, to the response. 
We follow the notation already shown in Figs. 1-3. 

Tables I and II give the results for the Spp channel. In Table I we compare all direct 
and exchange contributions from self-energy insertions. To avoid divergencies, diagonal self- 



energy insertions are evaluated up to infinite order |24|] . To do this, an average over the hole 
momentum of the bubble where the self-energy is attached should be done. This procedure 
is outlined in Appendix A. From the table it is clear that, although small in general, the 
exchange diagrams can amount to a non-negligible fraction of the direct ones, especially 
at energies around and below the quasiparticle peak. This is also visualized in Fig. 4, 
where the structure function including only the direct self-energy diagrams (long-dashed 
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line) is compared to that containing, in addition, the exchange ones (full line) for both the 
longitudinal (upper part) and the transverse (lower part) channels. The short-dashed line is 
the free structure function calculated with an effective mass of value of m*/m = 0.85. As has 



been observed before [|16| ,j32}-|34|l , the dressing of the nucleon lines by self-energy insertions 
smears out the structure function, moving strength out from the quasiparticle peak to the 
high missing energy region. 

Notice that in the diagrams shown in Figs. 1-3 we have not explicitly included the first 
order self-energy insertion on the fermion propagators. Instead of this, we have preferred 
to use a Lindhard function calculated with an effective mass m* in the nucleon propagators 
and, therefore, it already contains in an average way the effects of those self-energy insertions 
propagated to all orders. Actually, the use of m* is equivalent to having a real and energy 
independent self-energy parameterized by a function quadratic in the momentum of the 
nucleon. These mean field single-particle states define the HF ground state and the basis in 
which the perturbation theory has been constructed. 

The second order self-energy diagrams SE1D,SE3D,SE1E,SE3E are responsible for the 
appearance of an imaginary part which yields a width to the nucleon lines, which in turn is 
the responsible for the observed spreading in the nuclear response. 

In Table II we analyse the RPA-type correlations showing explicit results for the first and 
second order contributions. This is done for the transverse channel as RPA-type correlations 
are zero for the longitudinal one due to the election of our interaction with no / or /' 
Landau- Migdal terms. In the last column we also present the results for RPA correlations 
when exchange contributions are included up to infinite order, following the scheme of Ref. 



3T| . From Table II is clear that exchange terms of RPA type are very important. Their 
size is comparable (even bigger) to other direct diagrammatic contributions to the structure 
function and, therefore, they should not be neglected. Given the magnitude of these RPA 
exchange terms, the differences between the next-to-last and last columns of Table II also 
suggest that it is important to sum them up to infinite order as was done in Ref. The 
effect of the RPA diagrams in the transverse structure function is also displayed in Fig. 
5, where the full antisymmetric RPA series (full line) is compared to the direct ring series 
(long-dashed line). We observe that the transfer of strength from the low to the high energy 
region typical of the polarization (ring) diagrams is partly restored by the incorporation of 
the exchange diagrams. 

In Table III, we study the Sqp channel. In this kind of graphs the external operator 
creates (or destroys) a particle-hole pair and scatters a particle (or a hole). We have evaluated 
the case where the external operator scatters a particle. The case where it scatters a hole 
is negligible as can be found in Ref. pCfl . As mentioned above, due to the action of the 
antisymmetrization operator a graph where the external operator is attached to different 
bubbles has to be considered. This graph is SQPD' in Fig. 2 and has some influence as can 
be seen from Table III. The other exchange graphs are negligible. 

In Table IV the results for the Sqq channel are shown. As for the Sqp channel, exchange 
graphs are very small. The importance of the ground state correlation diagrams is clearly 
seen in Fig. 6, where the full lines represent the addition of Sqp + Sqq to the long-dashed 
lines, which contain the Spp contributions. In the longitudinal channel (upper part), only 
the self-energy terms contribute to Spp, while the transverse channel (lower part) contains, 
in addition, the RPA-type correlations. In view of these results it is clear that incorporating 
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the Sqp and Sqq channels is necessary in any perturbative calculation of the nuclear matter 
response as observed in Ref. [P0| . Having established in this work the smallness of the cor- 
responding exchange terms is particulary interesting, since the calculation of these channels 
can be restricted to the direct graphs thus avoiding a great deal of numerical computations. 

In addition our findings also support the idea that the use of an effective g' to account for 
exchange terms in the nuclear response is not appropriate. This point was already raised in 
Ref. PH|, where we built a prescription to calculate the full antisymmetric RPA series of the 



nuclear response. In that work, we showed that the use of a standard average prescription 
for g' was not able to reproduce the RPA antisymmetric response, especially for intermediate 
values of the momentum transfer. Using the average g' for calculating the other types of 
correlations would not be appropriate either because we have shown that they are basically 
dominated by the direct contributions. This is visualized in Table V where, changing the 
g' parameter to an effective value of g' = 0.5 does not, in the first place, reproduce the 
antisymmetric RPA response contained implicitly in the first column of the transverse part 
(a point already raised in Ref. |31|]) but also induces non negligible modifications in the other 
contributions, especially those related to ground state correlations (compare the next-to-last 
and last columns). 

From our study we conclude that the nuclear response is basically dominated by the 
direct diagrams and the exchange ones only need to be considered (and to all orders) for 
the RPA type correlations. 



IV. CONCLUSIONS 

A projection method to extract the main contributions of the ERPA theory with the 
explicit inclusion of exchange terms has been developed. This work is a continuation of a 



previous one |30[ in which only direct terms were studied. Here we have tried to clarify the 
importance of exchange terms of the particle-hole interaction, by performing a quantitative 
analysis of their influence in the nuclear response. 

The projection method classifies the contributions to the nuclear structure function into 
three channels, called Spp, Sqp and Sqq with P and Q being projections operators defined 
in Sect. 1. This separation permits to study the effects of the different types of correlations. 
In this sense, Spp represents final state correlations, Sqq ground state correlations and Sqp 
the interference between them. Through the analysis of our results we can conclude that all 
types of correlations are important and should be considered when one studies the nuclear 
response. 

After this statement, the the problem is the big number of graphs which should be eval- 
uated when the exchange part of the nucleon-nucleon interaction is retained. Before the 
numerical calculation we can see no reason to neglect any contribution. Our calculations 
show that for final state correlations, i.e. self-energy insertions and mainly RPA-type corre- 
lations, the exchange graphs are relevant in agreement with Refs. |37| , |38|| . On the other hand, 
they can be neglected for ground state correlations. Also, within the energy-momentum re- 
gion under consideration, exchange terms can not be parameterized by a redefinition of the 
g' parameter. 

It is also important to stress that the interaction employed and in particular the value 
for g', comes from parameterizations of processes at a lower energy-momentum region than 
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the ones considered here. Those values do not necessarily hold for us. Also from Table V, we 
see that a small change in one parameter can produce a noticeable change in the structure 
function. In any case, our scope was not a search of the optimal paremeters that produce a 
good agreement with the experimental data but a careful analysis of the exchange diagrams 
and our conclusions should remain valid in a wide variety of situations. 

In summary, our study shows that the nuclear response is basically dominated by the 
direct diagrams, the most relevant being those of Figs. 1-3, and the exchange contributions 
only need to be considered (and to all orders) for the RPA type correlations, which can be 
evaluated using the prescription of Ref . |HJ] . 
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APPENDIX A 



In this appendix we show explicit expressions for exchange contributions to the self- 
energy insertions Direct contributions can be found in Ref. Exchange self-energy in- 



sertions from eqs. (|24j) and ( Pq) (which contribute to graphs SE1E and SE3E of Fig. 1, 
respectively), are given by 

(z PQP (Q,v,h)r rtexch - = -j^- 4 (J^ fll^L J d ^ k Jm e(\h + Q-k\-i) 

6{l-\h + Q-k- k'\)9{\h + Q-k'\- l)Tl(k)Tl(k') 
(3g' 2 - (2(k.k') 2 - l)h' 2 + 2g'~h!) 

v-{Q 2 /2 + h.Q-k.k')+ir] ^ 

and 

(Re^ PRP (Q,^h)r^- h - = j^(J^ c fH^l j d * k JSk'9{l-\h + Q-k\) 

6{\h + Q-k-k'\- 1)0(1 -\h + Q- k'\)Tl(k)Tl(k') 
(3g' 2 - (2(k.k') 2 - l)h' 2 + 2g'h!) 

u-{Q 2 /2 + h.Q + k.k r ) ' ^ 

We have used dimensionless quantities Q = q/k F and v = Tiuj/2ef] k F and Bp being the 
Fermi momentum and energy, respectively. 

In order to simplify the calculation, it is a good approximation to eliminate the depen- 
dence on the hole momentum, by an average procedure (see Ref. |30H), as follows 

£ P( ^) p (Q, v) = J- [ d 3 h £ Pc ^) p (Q, u, h) (34) 

3 71 J 
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APPENDIX B 



In this appendix we show explicit expressions for exchange contributions to the structure 



function. Direct contributions can be found in Ref. |30 . 

Let us first consider the Spp channel. Graph SE2E of Fig. 1 (see Eq. fl2l"D ) is given by, 

(S sM ) L , T = (JL )> j d % j d 3 k I** OV tLnsE2E 

6(1 - \h\)0(\h + Q\ - 1)6(1 - \h + k'\)6(\h + k + k'\ - 1) 
6(1- \h + k\)9(\h + k + Q\ -1)(— Im)[ ' 



J u-(Q 2 /2 + Q.h)+ir] 

(35) 



v - (Q 2 /2 + Q.h + k.k') +irj v- (Q 2 /2 + Q.(h + k)) + irj 
Graph SE2E' of Fig. 1 (see Eq. (0)): 



(W(0,,)) W = ~ ( /a- %,n- 

0(1 - \h\)6(\h + Q\ - 1)6(1 - \h + k\)6(l -\h + k + k' + Q\) 
6(\h + k' + Q\)6(\h + k + Q\-l)(— Im)[ ' 



n u- (Q 2 /2 + Q.h) +ir] 
1 1 



, (36) 



v - (Q 2 /2 + Q.h - k.k') +ir] v- (Q 2 /2 + Q.(h + k)) + it] 

where 

OV L , S E2E = OV L , S E2E> = Tl(k) Tl(k') {3~g' 2 - (2(k.k>) 2 - l)h' 2 + 2~g'h'} , (37) 

and 

OVt,SE2E = OVt,SE2E' = 

= ) 2 rl(k) Tl(k'){A[h.(h + k) - (Q.h)(Q.(h + k)] 

[3g' 2 - (2(k.k') 2 - l)h' 2 + 2g'~ti] + 

(-3/i s 2 + fx v 2 )[Q 2 ~g' 2 + ((Q.k'f - 2(k.k')(Q.k)(Q.k'))h' 2 + 
(Q.k') 2 ~g'(k)h'(k') + (Q 2 - 2(Q.k) 2 )g'(k')~h'(k)} - 

(-3fi s + Hv)h' 2 -^ (k.k')[(k.Q)k' .(2h + Jfe) - (k'.Q)k.(2h + k)]} . (38) 

14 



First order exchange contribution to the RPA-type correlation (see Eq. (^) and graph 
RPAIE of Fig. 1): 

0(1 - fc)0(|/i + Q| - 1)0(1 - \h + k\)e(\h + k + Q\ - 1) 

(—- Im)U - - ) 

y 7T nK 2u- (Q 2 + 2h.Q) +ir) 2v + (Q 2 + 2h.Q) ' 

' ' )}, (39) 



2v - (Q 2 + 2(h + k).Q) +i7] 2v + (Q 2 + 2(/i + k).Q) 



where 



rl(k)(3g' + ti) (40) 

and 

OViTiRPAiE = (^ fTlmW + h') 4 + k) - (Q.h)Q.(h + k)]/Q 2 ] + 

(-3^ 2 + ^ > a )Q 2 (fif , + ^(fc.Q) 2 } . (41) 
Graph RPA2E of Fig. 1 (see Eq. (H)): 

| d 3 /i | d 3 A; CV (Li T),i?PA2£;0(l - h)6(\h + Q\ - 1)9(1 - \h + fc|) 
0(|/i + fc + Q|-l)( 1 ' 



2u - (Q 2 + 2h.Q) + ir] 2u + (Q 2 + 2h.Q) 
1 1 



2v - (Q 2 + 2(h + k).Q) +irj 2v + (Q 2 + 2(h + fe).Q) 



)}, (42) 



with 



£(Q, !/) = y rf 3 p + Q/2| - 1) 0(1 - \p - Q/2\) 

1 ' (43) 



2v-Q 2 -2Q.p + ir] 2v + Q 2 + 2Q.p 



and 
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OV {T ),rpA2E = Ti(Q)Tl(k)Q 2 Ag'^(g' + h'(k.Q) 2 ) . (44) 

(Note that for the present interaction the longitudinal contribution is zero). 
Graph RPA2E' of Fig. 1 (see Eq. ([23])): 



(S RPA 2E'(Q,v))l,t = ^ ( \*K j*k fM OV {l , m , 



9(1 - h)9(\h + Q\- 1)9(1 -\h + k\) 

6(\h + k + Q\ - 1)6(1 - \h + k + k'\)6(\h + k + k' + Q\ - 1) 

f_i T m )U - - ) 

K n nK 2u - (Q 2 + 2H.Q) +ir] 2v + (Q 2 + 2h.Q) ' 



2v - (Q 2 + 2(h + k).Q) +irj 2v + (Q 2 + 2(h + k).Q) 



2v - (Q 2 + 2(h + k + k').Q) +ir] 2v + (Q 2 + 2{h + k + k').Q) 



where 



and 



OV^rpme' = K(k) Ti(k')10[9(~g') 2 + h' 2 + 6g'h'} (46) 



ovpo.iu*^ = (^ (fc') {40[9(^) 2 + h' 2 + Q~g'~h'} 

[(h.(h + k + k')- (Q.h)Q.(h + k + k')]/Q 2 ] + 
9 J^l {Q 2 [(~g' 2 ) 2 + h' 2 (2(k.k>) 2 - 1)] + 

h' 2 [(Q.k) 2 + (Q.J?) 2 - 2(Q.k)(Q.J?)(k.J?)] + g' 2 h'[(Q.k) 2 + (Q.k 1 ) 2 }} . 

(47) 

Going now to the Sqp and S^Q-channel (see Eqs. ( p6f) and (|27|) ); we have for graph 
SQPE of Fig. 2 : 

(Ssqpe(Q,v))l,t = ^ ( / A / A / <ft< V (1 , W , 

0(1 - |/i|)0(|/i - fe| - 1)^(1^ - fe + Q\ - l)6(\h -k' + Q\-l) 

16 



0(\h + Q\- 1)9(1 -\h-k-k' + Q\)-rrnJ 



k.{k' - Q) 



( — Im)[ 



k v- (Q 2 /2 + Q-h - k.k') + irj v- (Q 2 /2 + Q.h)) + irj J ' 

(48) 



Graph SQPE' of Fig. 2: 



9(1 - \h\)9(\h -k\- l)9(\h -k' + Q\- l)9(\h + Q\ - 1) 

9(\h - k'\ - 1)9(1 -\h-k- k'\)-^— 

k.k 

^~ Im ^u- (Q 2 /2 + Q.h - k'.{k + Q)) + i v v - (Q 2 /2 + Q.h)) + i v ^ ' 

(49) 



where 



OV L , S QPE = OV L , S E2E , (50) 

Wt,sqpe = OVt,se2e , (51) 

OV t ,sqpe> = fTl(k) Tl(k'){4[h.(h + k') - (Q.h)(Q.(h + fc')] 

[3~g 12 - (2(k.k>) 2 - l)h' 2 + 2g'~ti] + 

(-3fx s 2 + fi v 2 )[Q 2 g' 2 + ((Q.k'f - 2(k.k')(Q.k)(Q.k'))h' 2 + 
(Q.k') 2 ~g'(k)h'(k') + (Q 2 - 2(Q.k) 2 )g'(k')h'(k)] - 

(-3fi s + fi v )h' 2 ^ (k.k')[(k.Q)k'.(2h + k') - (k'.Q)k.(2k + fc')]} . 

(52) 

Graph SQPZ/ of Fig. 2: 
(Sqpd'(Q,v))l,t = /«C* OV (I , nseFD , 

17 



6(1 - \h\)6(\h + Q\- l)6(\h - k\ - l)6(\h' + k + Q\-l) 
m + H - mi - m ^L film) 



v - (Q 2 /2 + Q.h) +i V u-(k 2 + Q 2 /2 + k.(h' -h) + Q.(k + h')) + irj ' 

(53) 

Graph SQQ3D' of Fig. 3: 

(wq,,» w - ^ (S-c ) 2 I3 J d3 '< l su l Sk 

6(1 - \h\)6(\h -k\- 1)6(1 - \h' + Q\)6(\h - k + Q\ - 1) 
6(\h' + k\ - 1)6(1 - \h'\) 



k 2 - k.h + k.h' k 2 + k.(h' -h) + Q.(h -h'-k) 

-1 1 
( V /m) v - (k 2 + Q 2 /2 + k.(h' -h) + Q.(h-k))+iri ' (M) 



In Eqs. (|53|)-(|54]), we have used, 

k' = k-Q 

and 

OV l ,sqpd> = OV L , S QQ3D> = rl(k) T 2 n (k') 5{Tg' 2 + (k.k'fh 12 + 2g'h'} , (55) 

OV T ,s Q pw = (^ ?Tl(k) Tl(k') 

{20[h.(h' + k) - (Q.h)(Q.(h' + k)}[3g' 2 + (k.k'fh' 2 + 2~g'~ti] + 

(3fi 2 + 2fi, 2 )[AQ 2 g' 2 + (Qx(kx k')) 2 h' 2 + 

(Q 2 + (Q.k') 2 )~g'(k)h'(k') + (Q 2 + (Q.k) 2 Yg\k')h'(k)] - 

(3fi s + 2fi v )h' 2 ^- (k.k')[(k.Q)(k'.(h + h' + k)) - (k'.Q)(k.(h +h' + k'))]} , 
Q 

(56) 

and 



OV T , SQQ3D , = ?Tl(k) Tl(k') 



{20[h'.(h - fe) - (Q.h')(Q.(h - k)}[3g' 2 + (k.k'fti 2 + 2g'h'\ + 

(3/i s 2 + 2fi v 2 )[4Q 2 g' 2 + (Q x (k x fc')) 2 /^ 2 + 

(Q 2 + (Q.P) 2 )s'(A;)/>'(A;') + (Q 2 + (Q.k) 2 )~g\k')h\k)} - 

(3/i s + 2fi v )h' 2 ^ (kJ?)[(k.Q)(k'.(h + h' - fc)) - (k'.Q)(k.(h + h' - fc'))]} . 

Q 

(57) 



Graph SQQ1E of Fig. 3: 

(S sqqMKt = ( ) 2 ^ /** 

OV (L> r),5ggi^(|p| - 1)0(1 -\p + k\)6(\p + Q\-1) 
6(1 - \h + k'\)9{\p + k + k'\ - 1) 



(fc.fc') 2 vr 7 V - (Q 2 /2 + Q.p + fc'.fc) + ir] 
Finally, graph SQQ2E of Fig. 3: 

0V (L , T ),5QQ2^(1 " - k\ ~ l)6(\h ~ k'\ - 1) 

6(\-\h-Q\)e(\-\h-k-k'\) 

1 Ira) [ - ^ 2/0 , ^ ,,,, , - ] , (59) 



where 



(k.k'f v 7T 7L z/ - (Q 2 /2 + Q.h - fe'.Jfe) + ir? J ' 

OVl,sqqie = OVl,sqq2e = OVl,se2e , (60) 

ov T , SQQ1E = ) 2 r 2 (fc) r 2 (A;0{[ P 2 - (Q. P ) 2 ] + qV 2 + /i, 2 ))/4} 

[3#' 2 - (2(fe.fc^) 2 - l)~ti 2 + 2~g'~h'] , (61) 
19 



and 



OV T , SQQ2E = ) 2 Tl(k) rl(k>){[h> - {Q.hf\ + + ^ 2 ))/4} 



>~ J2 {2(k.kf -l)h' 2 + 2g'h'] . (62) 
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TABLES 



TABLE I. Free and Self-energy contributions to the Longitudinal and Transverse struc- 
ture function (multiplied by 10 5 ). All results are for nuclear matter at momentum transfer 
q = 410MeV/c in units of MeV -1 fm -3 . The first column represents the energy transfer in 
MeV. Column Lind. represents the free structure function. Columns SE13 give the direct (D) and 
exchange (E) contribution to the diagonal part of the self-energy up to infinite order. Their first 
contributions are the graphs SEId, SE3d, SE1 e , and SE3 E of Fig. 1. Columns SE2 E , SE2e 
and SE2 E i are the non-diagonal self-energy contributions to the structure function as shown in 
Fig. 1. The last column is the sum of all these contributions given by Eq. (pi]). 



long. 


huj 


Lind. 


SE13 D 


SE13 E 


SE2 D 


SE2 E 


SE2 E > 


gLind.+SE. 


50. 


38.729 


-1.269 


0.715 


-1.570 


0.185 


0.071 


36.862 


100. 


46.106 


-4.390 


0.599 


-0.559 


0.374 


0.159 


42.289 


150. 


41.357 


-5.949 


0.397 


1.085 


-0.096 


-0.050 


36.743 


200. 


24.481 


-3.250 


0.111 


2.219 


-0.045 


-0.027 


23.490 


250. 


0.000 


2.971 


-0.127 


0.521 


-0.004 


-0.005 


3.357 


trans. 


50. 


60.380 


-1.938 


1.112 


-3.205 


0.037 


0.181 


56.566 


100. 


72.299 


-6.943 


0.946 


-0.604 


0.119 


0.004 


65.821 


150. 


64.610 


-9.354 


0.621 


2.575 


-0.071 


-0.232 


58.149 


200. 


37.739 


-4.948 


0.164 


5.337 


-0.082 


-0.357 


37.853 


250. 


0.000 


4.587 


-0.195 


0.726 


-0.006 


-0.079 


5.032 



TABLE II. RPA-type contributions to the Transverse structure function (multiplied by 10 5 ) 
in units of MeV -1 fm~ 3 for nuclear matter at momentum transfer q = 410 MeV/c. Columns 
RPAId (RPAl E ) and RPA2d (RPA2^ e+E i^) are the first and second order direct (exchange)-part 
to the RPA response, respectively. The notation is the same as in Fig. 1. Note that in that 
figure only forward-going contributions are shown while the present results contain both forward 
and backward-going contributions. Column RPA12 E > +E is the sum of all first and second order 
contributions (given by Eqs. (p2^) and (^)). Finally, column RPA an t. is the result for a full 



antisymmetric RPA using the formalism given in Ref. [31| 













RPA12 D+E 


RPAant. 


RPA1 D 


RPA1 E 


RPA2 D 


RPA2( E+E ,) 




50. 


-30.233 


8.844 


9.028 


-2.912 


-15.273 


-19.527 


100. 


-13.801 


4.209 


-3.098 


0.914 


-11.775 


-9.988 


150. 


5.122 


-1.132 


-3.441 


1.204 


1.753 


3.307 


200. 


10.354 


-3.206 


1.468 


-0.466 


8.150 


7.186 


250. 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 
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TABLE III. Longitudinal and Transverse Sop-type structure function (multiplied by 10 5 ) in 
units of MeV -1 fm -3 for nuclear matter at momentum transfer q = 410 MeV/c. The notation 
SQPd to SQPd' is the same as in Fig. 2. Column SQP an t is the sum of all contributions. 



long. 


hu 


SQP D 


SQP E 


SQP E > 


SQP D > 


SQPant 


50. 


4.950 


-0.003 


-0.077 


1.101 


5.971 


100. 


4.136 


-0.016 


-0.110 


1.142 


5.152 


150. 


3.116 


-0.055 


-0.075 


0.589 


3.575 


200. 


1.598 


-0.037 


-0.007 


0.234 


1.789 


250. 


-1.284 


0.004 


0.022 


-0.058 


-1.316 


trans. 


50. 


7.469 


-0.604 


-0.101 


1.324 


8.089 


100. 


5.412 


-0.672 


-0.144 


1.374 


5.969 


150. 


3.656 


-0.436 


-0.098 


0.928 


4.049 


200. 


0.906 


-0.001 


-0.009 


0.491 


1.387 


250. 


-3.652 


0.341 


0.029 


-0.073 


-3.356 



TABLE IV. Longitudinal and Transverse Sgg-type structure function. The notation SQQId 
to SQQ3d is the same as in Fig. 3. SQQ^Etot. is the sum of SQQ3e, SQQ3e> and SQQ3d' from 
the same figure. Column SQQ an t is the sum of all contributions. 



long. 


TlUJ 


SQQl D 


SQQl E 


SQQ2 D 


SQQ2 E 


SQQ3 D 


SQQ?>Etot. 


SQQant 


50. 


4.692 


-0.042 


2.455 


-0.074 


-0.552 


0.041 


6.520 


100. 


8.535 


-0.327 


1.428 


-0.046 


-0.629 


0.042 


9.003 


150. 


9.807 


-0.344 


0.661 


-0.019 


-0.299 


0.029 


9.834 


200. 


10.426 


-0.278 


0.310 


-0.008 


-0.220 


0.020 


10.249 


250. 


10.342 


-0.167 


0.047 


-0.003 


-0.113 


0.011 


10.118 


trans. 


50. 


6.064 


-0.051 


4.413 


-0.062 


-1.614 


0.012 


8.762 


100. 


12.673 


-0.307 


2.493 


-0.027 


-2.230 


0.017 


12.617 


150. 


15.153 


-0.313 


1.179 


-0.007 


-1.848 


0.014 


14.178 


200. 


17.352 


-0.263 


0.559 


-0.003 


-0.467 


0.012 


17.188 


250. 


16.551 


-0.167 


0.090 


-0.001 


-0.271 


0.006 


16.208 
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TABLE V. Comparison between our results and the corresponding direct values calculated 
with a modified value of the g'-Landau Migdal parameter {g' = 0.5) to partially reproduce the 
exchange contributions. Our results are given by SPP an t, SQP an t and SQQ an t where a value 
g' = 0.7 was employed with the explicit inclusion of exchange graphs. 



long. 


TlUJ 


SPPant 


SPPdir,g'=.5 


SQPant 


SQPdir,g'=.5 


SQQant 


SQQdir,g'=.5 


50. 


36.862 


36.409 


5.971 


5.034 


6.520 


2.529 


100. 


42.289 


43.256 


5.152 


5.442 


9.003 


4.995 


150. 


36.743 


36.430 


3.575 


4.312 


9.834 


5.971 


200. 


23.490 


19.059 


1.789 


0.966 


10.249 


6.291 


250. 


3.357 


2.393 


-1.316 


-1.336 


10.118 


6.015 


trans. 


TlUJ 


SPPant 


S PPdir,g'=.5 


SQPant 


SQPdir,g'=.5 


SQQant 


SQQdir,g'=.5 


50. 


37.039 


47.428 


7.109 


5.533 


8.762 


5.049 


100. 


55.833 


63.337 


5.362 


5.719 


12.617 


9.092 


150. 


61.456 


59.066 


4.049 


3.864 


14.178 


10.814 


200. 


45.039 


33.882 


1.727 


1.392 


17.188 


11.586 


250. 


5.032 


3.719 


-3.356 


-1.883 


16.208 


10.800 
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FIGURES 



FIG. 1. Goldstone diagrams stemming from Eqs. (^Tj)-(^). In every diagram the wavy lines 
represent the external probe with energy momentum (q, oS). The dashed line is the residual interac- 
tion. For simplicity we show only forward-going contributions, where the incoming external probe 
creates a particle-hole pair. In the backward-going diagrams (not represented here) the probe can 
also destroy a particle-hole pair. 



FIG. 2. The same as Fig. 1 but for the Sgp-channel given by Eq. (|2q). The action of the 
external probe represents the interference between scattering (a particle or a hole) and creating 
(or destroying) a particle-hole pair. 



FIG. 3. The same as Fig. 1 but for the S^Q-channel given by Eq. (27). The action of the 
external probe is to create (or destroy) a particle-hole pair. 
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FIG. 4. Self-energy contributions to the longitudinal (upper part) and transverse (lower part) 
structure function of nuclear matter at momentum transfer q = 410 MeV/c. Short-dashed line: 
Lindhard function using an effective mass m*/m = 0.85. Long-dashed line: effect of adding the 
direct self-energy terms. Full line: effect of adding the direct and exchange self-energy terms. 
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FIG. 5. RPA contributions to the transverse structure function of nuclear matter at q = 410 
MeV/c. Short-dashed line: Lindhard function (m*/m = 0.85). Long-dashed line: direct ring 
diagrams. Full line: full RPA structure function including the exchange terms to all orders. 
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FIG. 6. Contribution of ground state correlations to the longitudinal (upper part) and trans- 
verse (lower part) structure function of nuclear matter at momentum transfer q = 410 MeV/c. 
Short-dashed line: Lindhard function (m*/m = 0.85). Long-dashed line: Spp structure function. 
Full line: inclusion of the ground state correlation diagrams to the Spp structure function. 
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